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ABSTRACT 

A recent paper states that ultraviolet backradiation from the solar transition region and upper chromosphere strongly affects the degree 
of ionization of minority stages at the top of the photosphere, i.e., in the temperature minimum of the one-dimensional static model 
atmospheres presented in that paper. We show that this claim is incompatible with observations and we demonstrate that the pertinent 
ionization balances are instead dominated by outward photospheric radiation, as in older static models. We then analyze the formation 
of He in the above model and show that it has significant backradiation across the opacity gap by which Ha differs from other strong 
scatttering lines. 

Key words. Sun: photosphere - Sun: chromosphere - Sun: UV radiation 



1. Introduction 

Static one-dimensional (ID) modeling of the solar atmosphere 
assumes hydrostatic equilibrium in plane-parallel-layer geome- 
try to deliver temperature and density stratifications along ver- 
tical columns that may serve as reference description to com- 
pute spectral continua and lines similar to actual solar radi- 
ation. E arly versions of such "standard" models were com - 
piled by Ide Jageii (19591) and iHeintze et all (1964a| 1964bt . 
Holweger's (19671) LTE best-fit to optical iron lines became 
the first choice in classical abund ance determination, in partic- 
ular as the HOLMUL update by iHolweger & Miilleij (119741)1 
This model describes only the photosphere and is close to the- 
oretical photosphere models based on the additional assump- 
tion of energy conservation in the form of radiative equi- 
libriu m plus mixing-length convection in the deepest layers 
(e.g.. lKurucJ 1974 1994<; lGustafsson et alJ 1975, 2008). Semi- 
empirical best-fit modeling of solar continua including ul- 
traviolet wavelengths sampling the chromosphere continued 
with the Bilderbera Continuum Atmosphere (Gingerich & de 
Jager 1968), the Harv ard-Smithsonian Referen ce Atmosphere 
(JGingerich et al.l 19 711). the VAL model grids of IVemazza et al.l 
(il973i 1976' 1981') fitting disk-center Lyman brightness bins, 
the update by Maltbv et al. ( 1986), the more re cent update con- 
centra ting on disk-center ultraviolet spectra by lAvrett & Loeseil 
(1200 81) and the sequence of similar refinements bv Fontenla 

et al. ( 19901 19911 19931 20021 20061 20071 20091)^ 

This paper is triggered by the conclusion of iFontenla et al.l 



(I2009i) that "in the solar chromosphere, the FUV and EUV ob- 
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served emissions in continuum and lines produced in the upper 
chromosphere and transition region irradiate the low chromo- 
sphere and have significant ejfects on the ionization." . In partic- 
ular, they claim that such backradiation strongly affects the den- 
sity of minority stages of ionization in and near the temperature 
minimum of their quiet-Sun model. We first show that this claim 
is incompatible with basic observations, and then use illustrative 
computations to demonstrate that, instead, photospheric irradia- 
tion from below is the dominant age nt in ultraviolet continu um 
formation in the quiet-Sun model of iFontenla et al.l (120091)1 as 
was the case for the older standard models. 

We continue by demonstrating that significant chromo- 
spheric backradiation occurs in Htf in such static models through 
two-level scattering (Sect. |6]), but argue that the assumption of 
columnar hydrostatic equilibrium is untenable for the chromo- 
sphere and transition region (Sect. |7). 



2. Observation 

In this section we briefly review the observational character of 
the chromosphere, as counterpart to the standard-model defini- 
tion as the domain between the model's temperature minimum 
and steep temperature rise towards the corona. Figure [Tlgives an 
overview of the appearance of the solar atmosphere in an area 
containing plage and network. We use it to demonstrate that the 
actual chromosphere may have enough opacity to supply signif- 
icant backradiation in Hff, but not in the 1700 A continuum. 

The G-band image and the magnetogram in the first column 
of Fig. [T] are photospheric. The first image displays the gran- 
ulation and the spatial distribution of small kilogauss magnetic 
concentrations qualitatively and without sign through their proxy 
as intergranular bright points (best seen when zooming in). The 
second image quantifies the actual surface distribution of the 
line-of-sight component in such strong-field concentrations, but 
at lower angular resolution. The correspondence is excellent. 
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Fig. 1. Co-spatial simultaneous images from the Dutch Open Telescope (DOT) at 'ftp://dotdb.strw. leidenuniv.nl' and from the Solar 
Dynamics Observatory (SDO) at jhttp : //Imsal . com/get_aia_data| First row: DOT images sampling the G band, the core of Can H with 
1.4 A passband, the center of Ho- with 0.25 A passband. Second row: SDO/HMI magnetogram, SDO/AIA 1700 A image, SDO/AIA He II 304 A 
image. They were all taken on September 27, 2010 within a few seconds around 13:54:30 UT, a moment of excellent seeing at the DOT (Fried 
parameter ro ~ 14 cm). The field of view (scales in arcsec, terrestrial North up) is near disk center, South of the largest sunspot of AR11109. 
It contains a ridge of plage in its upper half, quieter areas in its lower half. The greyscales have been optimized by clipping (all images) and 
logarithmic scaling (AIA images). We invite the reader to zoom in with a pdf viewer for better appreciation of the detail in the DOT images. DOT 
observer: E. Romashets. DOT speckle reconstruction: A. Sukhorukov. Alignment routines: R Siitterlin. 



The other four images allegedly sample the chromosphere 
(1700 A, Can H, Ha) and the transition region (HeII304A). 
However, only the Ha image truly shows the chromosphere. It 
consists of long slender fibrils that cover and obscure the pho- 
tosphere everywhere in the form of overlying fibrilar canopies. 
This is always the case near active regions, plage, and active 
network; only the very quietest intern etwork areas are not cov- 
ered by fibrils in Ha (e.g., Plate 3.1 of Brav & LoughheadI 19741; 
iRouppe van der Voort et al.l 20071; iRutten et al.l 2008ll. These 
ubiquitous Ha fibril canopies constitute the ring of pink emis- 
sion just off the limb around the Sun that Lockver (1868) named 
"chromosphere". They are also likely to irradiate the low-opacity 
domain underneath with Ha backradiation. Their geometry re- 
mains unknown; the fibrilar appearance in images as this one 
may represent lona slender cylindrical fluxtubes (e.a., Heinzel 
& Schmieder 1994 ), density variations that cause corrugations of 
the T - 1 Eddin gton-Barbier surface in an extended atmosphere 
( Schoolmani 1 972 1). or warps in sheets as proposed bv lJudge et al.l 
(201 ID 

The 1700 A continuum shows no fibril canopies covering in- 
ternetwork areas but the underlying domain of cool gas that is 
every few minutes ridden through b y successive hot shocks pro- 
ducing Ca ii Htv and K2V cell grains ("Rutte n & Uitenbroekl 19911 
[Carlsson & Stein 1997, cf. Rutten 2011). It is called "clapoti- 
sphere" here following iRuttenI (1995) and is badly repres ented 
by the hydrostatic standard models (C arlsson & SteinI 19951). The 
bright points marking magnetic concentrations and making up 
the network and plage are roughly co-spatial with the deeper G- 
band bright points. Their large brightness contrast represents a 
mixture of scattered hot-wall radiation from below and current 
heating dCarlsson et al.li2010t) . Numerical simulations suggest 



that, due to the latter heating, magnetic concentrations do con- 
tain temperature stratifications resembling the standard models 
(e.g., the FAL-C model of Fontenla et al. 1993 ), but with the tem- 
perature minimum shifted down to height li^ 150 km, still within 
th e photosphere, due to fl uxtube evacuation (first quartet in Fig. 9 
of iLeenaarts et al.l 20101). Ubiquitous shocks occur also in and 
near these concentrations, but also already deeper th an in the 
shoc ked internetwork clapotisphere (ISteiner et al.lil998i Rutten 
et al. 201 II). Slow hydrogen ionization/recombination balancing 
in the cool aftermath of the shocks causes large NLTE overpop- 
ulation of the lower level of Ha (.Leenaarts et alj.2007b . 



The Ca 11 H image is much like the 1700 A image, apart from 
difference in angular resolution. In the internetwork it shows a 
bright mesh pattern that is dominated by shock interference and 
changes rapidly. There are no chromospheric fibrils. One should 
rather expect that fibrils that are opaque in Ha are also opaque 
at the center of Can H, because the thicker parts of Ha fibrils 
are seen also in Ca II 8542 A dCauzzi et al.l 2009l>. and Ca 11 H is 
bound to exceed Ca II 8542 A in extinction. Indeed, in the VAL3- 
C model Ca 11 K3 fo rms higher than the core of Ha (Fig. 1 of 
IVernazza et al.l 19811). However, filtergrams such as this one are 
made with wide passbands (1.4 A for the DOT) so that the in- 
ternetwork scene is dominated by shock brightness in the inner 
line wings that swamps the signature of very dark fibrils at line 
center (iReardon et al.l 20091). Compared to the 1700 A image, 
there is an extra contribution in the form of a diffuse bright haze 
around network and plage that is likely the unresolved appear- 
ance of t he highly dyn amic features seen as Ca n H straws near 
the limb (lRuttenli2006L . Ca n H spicules-II outside the limb (De 



R.J. Rutten et al.: Chromospheric backradiation 



Pontieu et al. I2007bt . and R BEs (rapid blue excursions) on th e 
disk in the blue wing of Ha (iRouppe van der Voort et al.l 2009t. 

The He II 304 A image shows fibrilar patterns that are not 
identical to but also not unlike the Ho- canopies. It also shows 
similar bright grains near fibril feet as Ho- does. These are not 
cospatial with the bright grains in the G-band and Ca ii H im- 
ages. The overall similarity suggests that hot transition-region 
gas shares the field-guided topography of the Ha chromosphere, 
outlining similar magnetic connectivity patterns. Such patterns 
seem also to be m apped by Lya fibrils in VAULT images 
(IVourlidas et al.l 2010). Any fibril that is opaque in Ha must be 
very thick in Lya but, comparably to Can H, the cell-covering 
fibrils are very dark in Lya whereas other structures, likely 
Doppler-shifted out of line-center obscuration, cont ribute larger 
brightness in the profile-summed VAULT images (iKoza et al.l 
2009). 

The conclusion from this section is that significant chromo- 
spheric backradiation may be expected in Ha and in He II 304 A 
which show opaque fibrilar canopies that constitute the actual 
chromosphere. Ca ii H should show them too at line center, but 
the DOT image and all other H & K filtergrams have too wide a 
bandpass. In the 1700 A image the absence of any fibrilar signa- 
ture, bright or dark or as erasure of the internetwork pattern seen 
in Ca II H, implies that they are transparent at this wavelength, 
even for such a not-so-quiet area. Since the backradiation by an 
optically thin fibril scales with the fibril opacity, it is unlikely 
that deeper layers are affected by chromospheric backradiation 
in such ultraviolet continua. 

3. Formulation 



iFontenla et al.l (I 200 9M based their claim of important chromo- 
spheric backradiation on their Fig. 4 containing graphs of pop- 
ulation departure coefficients b for selected levels of Si i. Si ii, 
Mg I, Mgll, Fe I, and Fe n. Population departure coefficients mea- 
sure the ratio of the actual population density (particles cm"-') to 
the population computed assuming LTE. There are two formats 
differing in normalizati on. The first, cal l ed "M enzel" here, fol- 
lows the original use bv lMenzel & Cillig ( 19371) who normalized 
hydrogen populations to the free proton density. The general- 
ization is to use a partial Saha-Boltzmann eval uation to obtain 
normalization to the next ion stage (E g. (6.3) of lJefferiesI 19681, 
Eqs. (13)-(15) of lVernazza et al.l I98II). The ojherf ormat, called 
"Zwaan" here, follows IWiiben ga & Zwaarj (19721) and uses the 
elemental abu ndance as normahzation. The two conventions are 
(lRuttenll20 03l): 



L,Zwaan 



= rii/n'r-''' and b\ 



Menzel 



tii/n 



LTE 



ric/n'^' 



LTE' 



(1) 



where n, is the population of level /, nc is the total population of 
the next ion, and the superscript LTE implies a complete solution 
of the Boltzmann-Saha LTE partitioning equations for the given 
level and all stages of the element. Sometimes the ion ground 
state ric is used instead of «c assuming «£. x nc- When the atom 
or ion containing level / is predominantly ionized, so that «c ~ 
«™ because most particles of the species are in the next-ion 
stage, then 



.Menzel _ ,, /„LTE 



but when most of the element sits in level / itself then 

T^Menzel ^ „LTE /„ ^ i /;,Zwaan 



«c /nc ~ 1/Oc 



(2) 



(3) 
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Fig. 2. T emperature strati fi cation in some standard models. Solid: 
FA L-C of [ Fontenla et "all (19931) Dot-dashed: its AL-C7 update 
by lAvrett & Loesej (20081) Dashed: FCHHT-B of Fontenla et al] 
(20091) Dotted: radiative-equilibrium model from Kurucz ( 1979 1992ai 
1 992b|) , extended to larger height assuming constant temperature. 



The continuum has b^™^''^ = 1 and the ratio bj/bc ~ b\jbc is the 
same in both definitions. 

Equations (|2l) a n d (|3]) were given as Eqs. (17) and (18) 
by Vernazza et aTl ( I98Ih who used the Menzel definition. 
IFontenla et al. f ( 2009^ did so too (private communication from 
J. Fontenla). We show results for both definitions here because 
the Zwaan version gives a more intuitive display of actual popu- 
lation departures for majority stages. For example, the Menzel- 
coeflicient dips for n - \ o f H i and C i at the tem perature min- 
imum in Figs. 30 and 33 of Vern azza et al.l ( 1981!) and Figs. 13 
and 14 o f Fontenla et al. ( 2007) do not imply that the populations 
of these levels are out of LTE. 



4. Demonstration 

We select magnesium f or plotting departure coefficients as in 
Fig. 4 of IFontenla et al.. (2009) Magnesium is an important 
electron donor in the upper photosphere where Mg i has larger 
Rydberg populations than anv other element (Fi g. 15 of Carlsson 
et al. 1992). It has important bound-free edges, at 162 1. 5 A (vac- 
uum) from the ground state (3s^ 'S") and especially at 2512.4 A 
(air) from the first excited level i^sip ^P). Mgi and Mgii spec- 
trum formation are characteristic for other abundant metals with 
low ionization energy (silicon, iron, aluminum). Together, these 
supply the electrons for the photospheric H opacity in the vis- 
ible and infrared and they dominate the solar c ontinuous opac- 
ity from the near-ultraviolet to Lya (Fig. 36 of lVernazza et al.l 
1981,). 

Figure [2] di s plays a selection of standard models. FAL-C of 
IFontenla et al.l (1 9931) was an update of VAL 3-C of Vernazza 
et al. (19811) with a less steep upper photosphere from the inclu- 
sion of the ultraviolet line haze and a different transition region 
from the inclusion of a mbipolar diffusion . A more recent update 
is the AL-C7 model of lAvrett & Loeserl (2008) constr ucted in 
particular to reproduce th e ultraviolet spectrum atlases of lBrekka 
(1993) and ICurdt et al. (2001) These spectra are also repro- 
duced by model FCHHT-B of IFontenla et alj (2009), which has 
a higher-located temperature minimum to accommodate dark in- 
frared CO lines. 

The dotted curve in Fig. |2]is a radiative-equihbrium model 
fro m Kurucz ( 1979 1992a 19921J). that we extended in Cauzzi 
et al. (20091) to larger height assuming constant temperature. It is 
very similar to the FCHHT-B model up to /; = 800 km, but has no 
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Fig. 3. Upper panels: magnesium ionization fractions A'siagc/Motai in 
the Kurucz model (left) and tlie FCHHT-B model {right), with A'stage 
the summed populations per ionization stage, A'lotai the total of all 
stages. Solid: NLTE. Dashed: LTE. Lower panels: gradient compar- 
isons. The y-axis units specify the logarithmic electron density A'e and 
total hydrogen density A'h in particles m"\ The temperature curves 
show 20 log(r) - 57 to provide comparable gradients. 



chromosphere or transition region and so serves as comparison 
for the case of no backradiation whatsoever 

We h ave used the ID version of the RH code (e.g., 
lUitenbroek 2001) in a setup with simil ar content as the 
PANDORA setup of lAvrett & Loesed (2008t) although the solv- 
ing method differs much. RH us es a multi-level accelerated 
lambd a iteration method following iRvbicki & Hum mer (1991 
'1992) - hen ce the code's nam e. We used the 66-level model atom 
for Mgi of ICarlsso n et al.l (1992) and a 12-level model atom 
for Mgii from lUitenbroek (1997) Mgiii is represented by its 
ground state only. Other elements of which the more important 
transitions are explicitly evaluated in NLTE are H, Si, Al, and 
Fe. Partial redistribution is accounted for in Mg nh&kandLyq 
and/? . The violet and ultraviolet line haze (e.g.. lGreve & ZwaanI 
1198 01) causing a quasi-continuum is accounted for by sampling 
the 236 000 lines between A - 1000 and 4000 A in the extensive 
list of Kumcf] every 20 mA, using a two-level coherent scatter- 
ing approximation per wavelength in which the collisional tran- 
sition probability of each sampled line is estimated from the ra- 
diative one with the approximation of Ivan Regemorterl (1962) 
rather t han applying a c ommon scattering recipe for all lines as 
done bv lAvrett & Loeser ( 2008) How e ver, a t est with the much 
simpler line-haze recipe of lBruls et al.l (ll992l)1 showed only little 
differences with the results shown here. 

Figure [3] shows the resulting NLTE and LTE ionization frac- 
tions for the Kurucz and FCHHT-B models. In both magnesium 
is virtually once-ionized throughout the atmosphere (below the 
transition region at /i = 1900km for FCHHT-B). Let us first in- 
terpret the simple LTE curves in the upper Kurucz panel, which 
reflect the competing effects in the Saha equation of the out- 
ward decreases of the temperature T and electron density A^e 
shown in the lower Kurucz panel. Both decreases are steepest 
in the deepest layers, where A^e drops from 10"' to 10""^ of the 
hydrogen density A^h because hydrogen becomes neutral. The 
Ne/Nii » 10"^ limit represents the combined abundance of the 
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Fig. 4. Population departures of magnesium ground states for the 
Kurucz model (upper panels) and the FCHHT-B model (lower panels). 
Left: Zwaan definition. Right: Menzel definition. Solid: Mgi ground 
state. Dashed: Mg ii ground state. Dot-dashed: Mg iii ground state. 



electron donor elements (Fe, Mg, Si, Al) and controls the H" 
opacity. The steep T drop wins from the steep A^e drop so that 
the Mg I fraction increases and the Mg m fraction decreases with 
height. Around h - 500 km the A^e drop wins from the T drop so 
that the LTE Mg i fraction drops slightly. The outer exponential 
decay of this fraction follows the decreasing gas density (with 
constant Ng/Nu ~ lO""*) at constant temperature. The NLTE Mg i 
curve also shows near-exponential decay around /! = 500km, im- 
plying that the nonthermal contribution, from radiative instead of 
collisional ionization and recombination, is about constant with 
height. 

The Mg I and Mg iii curves in the FCHHT-B panel are closely 
the same as in the Kurucz panel up to the FCHHT-B temperature 
minimum at h = 800 km, just as the models are. Thus, up to this 
height magnesium ionization does not depend much on higher 
layers. 

Figure |4] shows the population departures b of the magne- 
sium ground states for the Kurucz and FCHHT-B models. Both 
b definitions are used (left and right columns). Mg iii cannot be 
shown in the Menzel version because we do not include Mg iv. 
Since most of the population per ionization stage resides in the 
ground state, these curves reflect the divergences between the 
LTE and NLTE ionization fractions in Fig. [3] 

The FCHHT-B/Menzel curves in the last panel resem ble the 
corresponding curves in Fig. 4 of iFontenla et al.l (2009), which 
for Si II and Fe ii dip yet deeper than for Mg ii. At first sight the 
two deep dips in this panel would suggest large depletion of both 
Mg I and Mg ii at the temperature minimum. However, for Mg ii 
this dip is the reverse of the Mg iii peak in the Zwaan panel fol- 
lowing Eq. (O. In the Zwaan panels the Mg ii curves remain near 
unity because even large overionization of Mg i and overpopu- 
lation of Mg III do not much affect the Mg ii population, being 
dominant anyhow (Fig.O. 

Near and above the temperature minimum, the Mg i and 
Mgiii curves in the first panel of Fig. |4] mimic the temperature 
structure of the model atmosphere, for Mg iii reversely. This im- 
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Fig. 5. The Planck function By (solid), mean radiation Jy (dotted), and 
source function 5,, (dashed) averaged over the first 100 A of the IVlg i 
edges with thresholds at /1 = 2512.4 and 1621.5 A. Left: Kurucz model. 
Right: FCHHT-B model. Formal representative temperatures are plot- 
ted for these quantities: the electron temperature for By, radiation tem- 
perature for Jy, and excitation temperature for 5,,. They are obtained by 
applying the inverse Planck function to By, Jy, S ,,, respectively. The dot- 
dashed distributions at the bottom are the intensity contribution func- 
tions jy exp(-T,,), with jy the emissivity and t,. the optical depth,near 
the thresholds wavelengths. The lefthand deeper-formed curve is for the 
2512 A edge. Each curve is scaled to its maximum value. 

plies that the actual degree of ionization does not sense the tem- 
perature but is constant with height, so that the departure behav- 
ior is set by the LTE Saha-Boltzmann temperature sensitivity. 
A constant degree of ionization suggests dominance of constant 
radiation in an optically thin environment. 

5. Explanation 

The behavior of fo^waan ^^j. jy^g j ^^^ Fig. |4]is well-understood, in 
particular sin ce the classic Ph D analysis of the comp arable iron 
spectmm by ' LitesI ( 19721) (cf. lAthav & LitesI 1972; review by 
[Rutten 1988). The deep dip results from the dominance of the 
radiation field J in the source function S of the Mg i ionization 
edges in the violet and ultraviolet, where the bar denotes opacity- 
weighted frequency averaging over the edge and J is the angle- 
averaged intensity over all directions. 

Figure |5] illustrates this cause by plotting 100 A- wide aver- 
ages of Sy, By and 7,, for the two major Mg i edges against height 
in the Kurucz and FCHHT-B models. The spectral averaging 
smooths the many superimposed line blends. These quantities 
are plotted as formal temperatures in order to combine them in 
single diagrams for the two wavelengths, removing the Planck 
function sensitivity to wavelength (which is exponential in the 
Wien approximation valid here). Figure |5] shows that across the 
ultraviolet S y tends to follow Jy at the heights where the radia- 
tion escapes (shown by the contribution functions at the bottom 
of the graphs), with Jy x By for the 25 12 A edge. 

Single-transition scattering with S -(l-s) J+e B is an excel- 
lent approximation for these edges which scatter strongly, with 
complete redistribution over the edge profile. The steep decrease 
of the collisional destruction probability e with height makes S 
uncouple from B already deep in the photosphere, at a thermal- 
ization depth much deeper than the radiation escape depth. The 



sensitivity of the A operator to the steepness of the 5v(t) gradi- 
ent results in Jy > By where they uncouple. This excess produces 
a superthermal edge source function in the upper photosphere, 
with S s» (bc/bi) B in the Wien limit for an edge between level 
/ and ionization limit c. Thus, the divergence between Sy and 
By for the edge at 1621.5 A in Fig. |5] translates into the diver- 
gences between be (Mgii) and bi (Mgi) in Fig.|4]that measure 
the source function departure from LTE for this edge. Similarly 
for bc/b2 in the edge at 2512.4A (which is a more important 
ionization channel and continuum provider). 

In a plane-parallel atmosphere Jy flattens out to a constant 
outer limit unless there is significant emissivity in higher layers. 
The ID models have that in their transition regions where the 
steeply increasing temperature gives larger s, causing renewed 
coupling to By. Below the transition region, the flat Jy may cut 
through whatever the model specifies as chromospheric temper- 
atures, so that the latter are mapped into the bc/bi ratio when Sy 
follows Jy. In the FCHHT-B panel of Fig. |5] the 25 12.4 A edge 
shows such behavior. The 1621.5 A edge has some Sy coupling 
to the higher chromospheric temperature and a yet smaller in- 
crease in Jy, which rises slightly above the flattened-out limit 
value for the Kurucz model. 

All this i s beautifully illustrated across the VAL3-C spectrum 
in Fig. 36 of Vernazza et al.l (iJWlhi and is explained at length in 
iRuttenl l 20031) 

There is no significant difference between the two models 
in Fig, |5] at the heights of formation of these bound-free Mg i 
continuum contributions. For both edges the formation is purely 
photospheric. Backradiation from the upper chromosphere plays 
no role. 



The conclusion from Figs. |2]-|5]is that the FCHHT-B model 
differs quantitatively from VAL3-C, but not qualitatively and not 
in its radiation physics. In particular, the continuum formation at 
these ultraviolet wavelengths is identical between the two mod- 
els used here since the diff'erences in 5 v at larger heights are not 
sampled by the photospheric contribution functions. 

Figure |6] adds FCHHT-B b curves for selected excited levels 
in Mg I and Mg n to the ground-state display in the upper panels 
of Fig. |4] The Mg i ground state has the largest NLTE deficit in 
the temperature minimum due to J > B overionization in the 
principal Mg i edges. The plateau of bi (index 0) excess at larger 
height maps J < B in these scattering edges. 

The higher levels in Mg i converge to the b^l value of the 
Mgii ground state, with multi-level crosstalk ("interlocking"). 
In the Mg i Rydberg domain this convergence is maintained by 
a collisional population replenishment flow driven by photon 
losses in intermediate Mg i lines, exemplified by the deep onset 
of the drop in the « = 10 curve. Such photon-loss replenishment 
fro m the next-ion pop ulation reservoir was called photon suction 
bv lBruls et al.l (19921) The Rydberg flow causes the Mgi emis- 
sion features near A - 12/im that were thought chromospheric 
bylZirin & Podd ( 1989) but actuallv are pho tospheric (Carlsson 
etal. 1992). 

The Mgii b peaks in the second panel of Fig. |6] reflect the 
effect of J > B pumping in strong ultraviolet lines, similarly as 
in Fell (see iRutten 1988). A particular example of the latter is 
the weak Fe n 3969.4 A li ne between Ca ii H and He, discovered 
as limb emission line by lEvershedl ( 1 9291) Its on-disk emission 
was initially attributed to chromospheric backradiation bv lLitesI 
( 1974), but is ac tually caused by p hotospheric pumping in Fe ii 
resonance lines dCram et al.l 1980) which ex plains its extraor- 
dinary spatial intensity variation at the limb (iRutten & Stencell 
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Fig. 6. Population departures of selected magnesium levels in the 
FCHHT-B model. The level numbers increase with excitation energy. 
Left: Mgi. Right: Mgii. Upper row: Zwaan definition. Lower row: 
Menzel definition. 



I198CI I. Similar upper-photosphere J > B pumping sets the b peak 
of the Mg m ground state in the first panel of Fig. |4] 

The lower panels of Figure |6] repeat the b displays for the 
Menzel definition. The lefthand panels are the same, but on the 
right the Menzel panel shows reversals of the Zwaan panel due 
to the normalization by the Mg m population. It makes the Mg ii 
curves appear similar to the Mg i curves. 

The lower panels of Figure [6] are similar to the correspond- 
ing Mg I and Mg ii panels in Fig. 4 of lFontenla et al.l (2009) The 
look-alike sequences of dips, deepest for the lowest levels, in 
their six Menzel plots made Fontenla et al. (2009) conclude: "All 
our results display overionization (i.e., ground-level departure 
from LTE coefficients smaller than unity) of all species around 
the temperature minimum. [. . . ] Also, near the temperature min- 
imum the lower levels have smaller departure coefficients than 
the upper levels, indicating that the overionization of neutrals 
is a result from FUV and/or EUV irradiation which primarily 
affects lower levels. This indicates that overionization is much 
more affected by the irradiation from the upper chromosphere 
in continuum and emission lines than by the photospheric ra- 
diation and absorption lines. Therefore, we stress that the con- 
sideration of a realistic upper chromosphere is essential to the 
determination of the densities near the temperature minimum of 
(1) neutral low first-ionization-potential (FIP) elements, and (2) 
singly ionized high FIP elements. Even a very sophisticated cal- 
culation of the effects of lower chromospheric absorption lines 
on the elemental ionization can produce unrealistic results if it 
does not include upper chromosphere irradiation." . This inter- 
pretation is incorrect. Downward irradiation does not aff'ect the 
degree of ionization of these minority species around the temper- 
ature minimum of the FCHHT-B model. We have demonstrated 
this for Mg here; the same applies to Fe, Si, Al. 

However, more positively, we wish to note that the model- 
ing of iRintenla et al. (2009) is not impaired by this misinter- 
pretation of its results, and that that does not detract value from 
their effort. In this era of giant computer programs it becomes 
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Fig. 7. Ho- formation for the Kurucz (left) and FCHHT-B (right) mod- 
els. Solid: Planck function By at line center. Dashed: angle-averaged ra- 
diation Jy at line center. Dotted: total source function 5 ,, at line center. 
The ticks mark the locations with line-center optical depth Ty = 3, 1, 0.3 
from left to right. 



non-trivial to diagnose what the output implies. Chromospheric 
physics is becoming a key area in solar research but does require 
understanding of complex non-equilibrium optically-thick spec- 
trum formation. Our tutorial elucidation of this particular issue 
may be instructive to newcomers to the field. 



6. Backradiation in Ha 

Figure Q shows formation parameters for H o; in the Kurucz and 
FCHHT-B models. It is similar to Fig. 8 in lCauzzi et aP (2009|U 
but here, prompted by the referee and in the vein of this paper, 
we add detailed analysis using the FCHHT-B model for demon- 
stration. We do not regard it a viable explanation of actual solar 
He formation (Sect. |7]i, but use it rather as a one-dimensional 
didactic "FCHHT-B star". 

The lefthand panel of Fig. |7] illustrates Ho- line formation 
in such a star without overlying chromosphere. At line center, 
Jv drops well below By. This behavior differs from the Jy > By 
excesses of the ultraviolet continua in Fig. |5]because the Planck- 
function sensitivity to temperature is smaller at longer wave- 
lengths, flattening the dBy/dh gradient, and because the substan- 
tial additional line opacity a'y flattens the dBy/dTy = dBy/Ka^ + 
al) dz] gradient yet more. The A operator then produces Jy < Sy 
tending towards the S v(0) = yfs^By surface value for a scattering 
isothermal atmosphere with constant Sy. 

In the Kurucz model the Hff core originates from the upper 
photosphere around h = 400 km as a deep, strongly scattering 
line with ^v ~ ^vo- The corresponding emergent profile is shown 
in the top panel of Fig. [8] It is well reproduced by applying the 
Eddington-Barbier approximation /,, ^ Sy{Ty-l). 

The righthand panel of Fig. [T] shows the formation of Ha 
in the FCHHT-B model. Its chromospheric high-temperature 
plateau supplies sufficient Hff opacity that at line center r,, = 1 is 
reached about a thousand km higher than in the Kurucz model. 
The line-center formation is again well described hy ly ^ S y{Ty = 
1) ~ -/vo(tv = 1), but the Eddington-Barbier approximation is less 
well applicable in the line wings which hav e a formation gap 
between photosphere and chromosphere (e.g.. lSchoolmanll972t 
iLeenaarts et al.l 2006l). This gap is detailed in Fig.|9]which shows 
Hff opacities and optical depth scales, with the opacities multi- 
plied by the scale height in the upper FCHHT-B photosphere (it 
doubles in the FCHHT-B chromosphere from a steep increase of 
the imposed non-gravitational acceleration). The line extinction 
has a deep dip in the temperature minimum. Correspondingly, 
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Fig. 8. Top panel: emergent Ho- intensity profiles at disk center com- 
puted from the FCHHT-B (dotted) and Kurucz (dashed) models, in 
comparison with the observed spatially-averaged profile (solid) taken 
from t he atlas observed by Brault & Testerman. calibrated by Meckel 
& Labs ( 19841), and posted bv lNeckell ( 19991) Lower panels: mean in- 
tensity Jy across Ho- at diff'erent heights in the Kurucz and FCHHT-B 
atmospheres, respectively. For the Kurucz panel the sampling heights 
are 26, 41, 157, 283, 358 and 995 km. For the FCHHT-B panel they are 
35, 150, 300, 519, 850, 1000, 1340 and 1701 km. The dotted profile in 
the bottom panel is the outward intensity at height 850 km, impinging 
on the FCHHT-B chromosphere. The dotted profile in the middle panel 
is the outward intensity at height 252 km in the Kurucz model where 
the line-center optical depth equals the optical thickness t,, = 3.5 of the 
FCHHT-B chromosphere. 



the optical depth buildup (dotted curves) levels out. At line cen- 
ter the FCHHT-B chromosphere has optical thickness Ty -3.5. 
At A/l = -0.44 A Ty = 0.3 (outer tick) is reached already in the 
chromosphere but Ty = 1 only a thousand km deeper in the pho- 
tosphere. At A/l = -0.74 A the FCHHT-B chromosphere is vir- 
tually transparent. 

At line center the line extinction much exceeds the continu- 
ous extinction even in the dip, so that the total source function 
equals the line source function at all heights. In the wings the 
line extinction drops below the continuum extinction. Where this 
happens the total source function drops to the Planck function 
(not shown). 

Figure [To] shows the departure coefficients ^2 and b^, (Zwaan 
definition) governing Ha in the FCHHT-B model. The curve di- 
vergence logib^) - log(/72) corresponds to the curve divergence 
logiSy) - logiBy) in the FCHHT-B panel of Fig.]?] Throughout 
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Fig. 9. Solid: Ha line extinction coefficients a[ per m against height for 
the FCHHT-B model, at line center and in the blue wing at A/1 = -0.44 
and -0.74 A, and the continuous opacity a^ in the adjacent continuum. 
The dot-dashed curve is the line-center result for LTE. The opacities 
are multiplied by a scale height of 10' m to make them comparable to 
the corresponding optical depth scales (dotted curx'es). The tick marks 
on the latter conespond to r,, = 3, 1, 0.3. Across the opacity dip due to 
the temperature minimum the optical depth is Ty = 3.5, 0.7, 0.04 at the 
three line wavelengths, respectively. 
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Fig. 10. Hff population departure coefficients against height for the 
FCHHT-B model. The tick marks correspond to Ty = 3, 1,0.3. Solid: 
lower-level departure coefficient 62- Dashed: upper-level departure co- 
efficient ^13. 



the near-isothermal FCHHT-B chromosphere the Ha opacity is 
nearly in LTE because the formation of Ly a is close to detailed 
balancing there. However, this is not the case in the deep tem- 
perature minimum and steep rise to the chromospheric tempera- 
ture. Here the A operator smooths the rapid temperature changes 
for the scattering-dominated Lya source function (not shown) 
which has 5v ~ /v ~ ^2 By. This Lya smoothing causes the high 
peak in b2 and considerable corresponding fill-in of the opacity 
dips in Fig. |9] illustrated for line center as difference with the 
much deeper dot-dashed LTE curve. In this manner Lya scat- 
tering compensates for low-temperature Ha opacity loss when 
cool features are narrow in spatial extent. The b2 peak in Fig.fTOl 
doubles if CRD is adopted for Lya, from additional smoothing 
through farther wing-photon travel. 

The corresponding emergent intensity profile in the top panel 
of Fig. [8]is similar to the Kurucz result despite the disparate line 
formation. Both are reasonably good approximations to the ob- 
served mean quiet-Sun profile. Both gain fit quality in the wings 
if the computed continuum intensity is rescaled to the observed 
value. Increase of the microturbulence in the Kurucz model from 
1.5 to about 10 km s"' produces a good core fit for this model as 
well. 
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The Jy curves in Fig. [T] show that the opacity gap in the 
FCHHT-B model contains much more Ho- radiation than at cor- 
responding heights in the Kurucz model. This Wa radiation fill- 
in causes a yet higher peak for b-^ in Fig. [10] It is detailed in 
the lower panels of Fig. [8] In the lower photosphere the Ha Jv 
profiles are the same for the two models, but above h ^ 250 km 
the Kurucz model gives absorption cores whereas the FCHHT- 
B model predicts slightly self-reversed cores across the opac- 
ity gap. This difference explains that for the Kurucz model Jy^ 
(shown by the line-center Sy curve in the lefthand panel of 
Fig. |7]l ends up higher than Jy at line center, whereas for the 
FCHHT-B model Sy =» 7v„ coincides with line-center Jy across 
the opacity gap in the righthand panel of Fig. |7] The large 
core width of the highest-formed FCHHT-B profiles comes from 
large thermal broadening (about 10 km s"') and yet larger mi- 
croturbulence (about 15 km s"'). 

Because the FCHHT-B chromosphere is near-isothermal the 
Ho- formation is comparable to the classic results for a finit e 
isothermal scattering atmosphere of lAvrett & Hummer! (I965h 
in particular their optically thick but effectively thin case in 
which the line source function does not reach thermalization. 
The main difference is that the FCHHT-B chromosphere is ir- 
radiated from below. The impinging intensity profile in the out- 
ward direction is shown by the dotted curve in the bottom panel 
of Fig. [8] This line is much shallower than the emergent pro- 
file from the Kurucz model, in agreement with the much higher 
Jy above h ^ 200 km. A test with a hotter chromosphere, i.e., 
thicker in Ha, gives more backradiation and a higher imping- 
ing profile with a bright self-reversed core. However, the deep 
emergent-intensity core remains about the same. 

Thus, the presence of an opaque chromosphere changes the 
illumination profile with respect to the profile that emerges from 
the same photosphere without overlying chromosphere. A math- 
ematic al explanatio n is given by the principle of invariance (page 
165 of IChandrasek har 1950 1 regarding addition or subtraction 
of a layer of arbitrary optical thickness to a semi-infinite plane- 
parallel atmosphere; creating a gap or pushing down the finite 
chromospheric atmosphere over the gap to where it meets simi- 
lar source conditions does not change the emergent radiation. A 
more physical explanation is that the radiation field in the gap 
builds up due to bac kradiation from the ov erlying chromosphere 
(called reflection bv lChandrasekhaiJ 1950ll. 

We now demonstrate that this backradiation is dominated by 
two-level scattering. The line source function for complete re- 
distribution can be written as: 



(1 



■ 77vi,) 7vo + Sv„By^{T) + rjy^By^(Ti) 



Jvo+<,By„(T) + Tly,By,(Ta) 

1 + 4„ + '?;„ 



(4) 



(5) 



where Jy^ is the profile-averaged mean intensity, Sy^ the thermal 
destruction probability, i.e., the fraction of line-photon extinc- 
tions by direct collisional deexcitation corrected for emissivity 
from spontaneous scattering, rjy^ the fractional probability of all 
indirect detour extinction, with detour meaning any multi-level 
path from the upper to the lower level not including the direct 
transition, corrected for stimulated detour emissivity, and ej,^ and 
TyJ,^ the corresponding ratios of such extinctions to the contribu- 
tion by two-level scattering. The formal detour excitation tem- 
perature Td is given by (guDui)l(giDu,) = sxp{livo/kTd) where 
£)„/ is the summed transition probability (per second per parti- 
cle in the upper level) of all detour upper-to-lower paths, D/„ for 
all detour lower-to-upper paths, and gi and gu are the statistical 
weights. 
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Fig. 11. Ha source function comparisons for the FCHHT-B model. 
Solid: Planck function and total source function at line center The lat- 



ter has ticks at 
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3, 1,0.3. Dotted: line source function 5' , dif- 



fering from the total source function only at the very top of the atmo- 
sphere. Dashed and dot-dashed: two-level approximation ( 1 - £,,„ ) y,,„ + 
BygBy^iT) and profile-averaged mean intensity Jyg. These coincide ev- 
erywhere. 
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The classic literature used Eq. |5] (e.g., iGebbie & Steinitz' 
1974, Sect. 8.1 of Jefferi es 1968, Eq. 12.11 of Mihalas 1970) 
following iThomasI (119571)1 who divided strong lines into "colli- 
sion type" with Syg > rjyg and "photoelectric type" with ^vo > ^vo' 
designating Ha as principal example of the second type. In par- 
ticular. Ha would gain most photons from Balmer ionization 
followed by a recombination path into » = 3 employ ing pho- 
tospheric Ba lmer photons (e.g.. lAthav & ThomasI 19611 see also 
ICramI 19851). The correspon ding schemati c sourc e function di- 
agram in Fi g. 3 o f Jefferie s^ & ThomasI (19 59h (reprinted in 
Fig. 12-9 ofiMiiialas 1970 and Fig. 11-11 of lMihalasI 1978') is 
qualitatively similar to the FCHHT-B Ha behavior in Fig.|7] with 
S 'y leveling out in the photosphere to become much higher than 
By in the temperature minimum. 

However, in the FCHHT-B case this behavior is due to 
the backscattering from the chromosphere, not from Balmer- 
continuum detour emission. If the latter were important it would 
raise the Ha source function also for the Kurucz model. The 
Balmer continuum scatters outward similarly to the 2512 A con- 
tinuum in Fig. |5] the same in both models up to h x 800 km, 
with about 1400 K radiation temperature excess in the tempera- 
ture minimum - but that region is transparent in Ha. 
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The scattering nature of the Ha source function is demon- 
strated in Fig. [TT]by comparing it with the mean radiation /vo 
and the two-level approximation. They are closely the same 
everywhere. The differences are magnified in Fig. [12] showing 
the fractional collision and detour contributions. The collisional 
contribution is negligible except in the deep photosphere where 
collisional excitations create most Ha photons (with Sy^ =s 1 and 
Jvf, ~ Bvo)- The latter resonance-scatter outward. Higher up, the 
detour contribution is m uch larger, confirming the riy^/sy^ evalu- 
ation of lAlet all (20041) but it still amounts to only a few percent 
of the line source function and raises the line-center intensity by 
only one percent of the continuum intensity. In the steep temper- 
ature rise of the FCHHT-B transition region the detour contribu- 
tion grows rapidly to nearly 100%, but at negligible Ha opacity 
and therefore of no importance to the emerging profile - and not 
due to photospheric irradiation. 

The conclusion from this section is that in the FCHHT-B 
model Ha is first and foremost a scattering line, with chromo- 
spheric backscattering boosting the line-center radiation across 
the opacity gap between photosphere and chromosphere. The 
Ha core emerges from the FCHHT-B chromosphere but most 
photons were created in the deep photosphere. For Ha the 
FCHHT-B chromosphere is primarily a scattering attenuator that 
builds up its own irradiation from below. The main source func- 
tion difference with other strong scattering lines such as the Ca ii 
resonance lines and infrared triplet is the effect of the opacity 
gap, not the amount of detour contribution. 



7. Discussion 

Static ID modeling of the continua from the solar atmosphere 
has reached such sophistication that even the modelers them- 
selves may misinterpret their results (Sect. |5)- Interpretation of 
chromospheric fine structure using strong lines, in particular Ha, 
has relied more on cloud modeling (see the excellent review by 
iTziotzioul 2007) than on static model-atmosphere interpretation 
as in Sect.|6] In the meantime, observations of the actual chromo- 
sphere (Sect. |2]i have reached such spatial and temporal resolu- 
tion that the worlds of the modelers and the observers now seem 
far apart. A third worldview is given by time-dependent numeri- 
cal MHD simulations of the chromosphere (e.g .^Schaffenberger 
et al., 2006.: .Leenaarts et al. 2007 2010, Martinez-Svkora et all 
1201 H . Obviously, the different views should come together It 
falls outside the scope of this paper to effectuate such synthesis, 
but we make a few salient points. 



Emergent sp ectrum fitting. The FCHHT-B model and the AL- 
C7 model of lAvrett & Loeser (2008.) fit the ultraviolet atlases 
about equally well, but differ appreciably in their stratifications 
(Fig. |2| and in their stratification physics. Both assume hydro- 
static equilibrium, but the AL-C7 model gains additional chro- 
mospheric support from imposed turbulent pressure inspired by 
observed line broadening, as in its VAL and FAL predecessors. 
The FCHHT-B model gains chromospheric support instead from 
imposed non-gravitational acceleration inspired by the Farley- 
Buneman instability (Fontenla et al. 2008 ). Thus, best-fit model 
determination from such data does not represent a unique data 
inversion. 

An extreme example of non-uniqueness is that the two dis- 
parate Ha formations in Fig. |7] with and without a chromo- 
sphere, can both fit the observed average Ha profile (top panel of 
Fig.[8]l. Figure [T]obviously suggests that chromospheric Ha for- 
mation as in the FCHHT-B model is more realistic than photo- 



spheric Ha formation as in the Kurucz model, but the predicted 
Ha profile cannot be used as discriminator. 

One-dimensional continuum modeling. The similarity of the 
magnetogram, 1700A and Can H panels in Fig. [1] suggests that 
ID modeling along vertical columns which differs between in- 
ternetwork, network, and plage may be a reasonable approxi- 
mation for ultraviolet continua from these features. The bright 
grains that constitute network and plage mark strong-field struc- 
tures that are roughly vertical. However, facular modeling to- 
wards the limb must then account for slanted viewing across 
such structures by cutting through different ID models along the 
line of sight, as in the classic Zurich fluxtube modeling of e.g., 
lBunteetalJ(.1993t)l 



Cloud modeling. The FCHHT-B chromosphere may be re- 
garded as an horizontally extended isothermal slab but not as a 
constant-property cloud since it is more like a finite atmosphere, 
with density stratification, outward Sy decay, and Eddington- 
Barbier-like Ha core formation. The source function is dom- 
inated by scattering (Fig. fTTT i and the slab is effectively thin 
(Fig.|9]i so that the source function amplitude depends on the slab 
opacity. Due to backscattering, the Ha profile illuminating the 
slab from below (dotted curve in the bottom panel of Fig. |8]l is 
shallower than for locations without an overlying chromospheric 
blanket (Kurucz prediction in the tot) panel of Fig.lSl cf. Bostanci 
& Al Erdogan 2010) and brightens for a thicker blanket. If actual 
fibrilar canopies share such properties than these upset classical 
cloud modeling. 

At the suggestion of the referee, we added in the middle 
panel of Fig.[8]the outward intensity profile at that optical depth 
in the Kurucz model which equals the optical thickness of the 
FCHHT-B chromosphere. It is similar to the outward intensity 
profile impinging the latter (bottom panel) because the line for- 
mation is much more similar for the two models on optical depth 
scales than on geometrical height scales. The near-equality of 
these profiles suggests a new recipe for ID Ha cloud model- 
ing, namely to use the outward intensity profile in a radiative- 
equilibrium photosphere at optical depth equal to the cloud's 
thickness as impinging background profile. 



Dynamics. The actual internetwork atmosphere above the pho- 
tosphere (in our view consisting of a cool clapotisphere, chro- 
mospheric fibril canopy, and a similar-morphology sheath-like 
transition region) is continuously shocked. In Ha ubiquitous 
internetwork s hocks are evident in Dopp ler timeslices (e.g., 
iRutten et "all 20081; ICauzzi et"aLl 20091). In magnetic con- 
centrations shocks are ubiquitous already in the photosphere. 
Shoc ks also p roduce the d ynamic fibrils jutting o ut from net- 
work (|Hansteen_etaL 2006, Dejo^euietdj2007a|). The elusive 
straws /spicules-II/RBEs near network are even more dynamic. 
This basic chromospheric dynamism upsets static chromosphere 
modeling. 



Departures from LIE. The FCHHT-B NLTE departures dis- 
cussed here reach two to three orders of magnitude at the 
FCHHT-B temperature minimum (Figs. |4] [TOl i while the lower 
level of Ha is nearly in LTE where the line forms (Fig. [TOt . The 
non-equilibrium simulation of Leenaarts et a l. (2007) predicts 
overpopulations of this level up to twelve orders of magnitude 
for cool post-shock phases, where the temperature may fall well 
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belo wthe temperature minima of the standard models ( Leenaarts 
et al. 201 11). Thus, dynamic behavior may dramatically increase 
the temperature variations and NLTE effects. 



Height of the transition region. One-dimensional models for 
different features typically differ primarily in t he height of their 
transition region. Simulations as the one of Leenaarts et alJ 
(2007i) suggest that the transition region above internetwork is 
kicked up and mass-loaded by field-guided flows that do not 
obey one-dimensional hydrostatic equilibrium. In this simula- 
tion the location of the transition region varies over h = 1000 - 
4000km, with rapid changes, and is generally lower above mag- 
netic concentrations. Such larger range than in the ID model 
grids may better reflect actual variations in fibril-canopy height 
or corrugation. 

Ha opacity gap. The Ha opacity gap and backscattering into 
it are properties of the FCHHT-B model. Do they also occur in 
the real Sun? The images in Fig.[T]roughly agree with FCHHT-B 
predictions: in the 1700 A image no fibrils are seen whereas they 
appear opaque in the Ha image. Also, when one samples the 
real Sun in Ha away from line center (e.g., in DOT movieo the 
photospheric granulation appears when the fibrils become trans- 
parent, without an intermediate clapotispheric scene as in the 
Can H and 1700 A images in Fig.[T] Thus, a similar opacity gap 
seems to exist under the actual fibrilar canopies making up the 
internetwork chromosphere. The radiation crossing it, boosted 
by backscattering, will suffer substantial spatial smoothing of 
the photospheric scene. Some of it may leak out through an ef- 
fectively thin fibril canopy, s catter aro u nd dar k fibrils in photon 
channeling as suggested bv I Al et aP (20041), or be seen from 
aside as in the backradiation explana tion of bright rims under 
filaments (e.g.. lKostik & Orlovali 197511 cf. lPanasencq 20101) . 

Ha detour brightening. In the FCHHT-B model Ha is an al- 
most pure scattering line (Fig. fTTt . The detour contribution 
T]vo(Bvo{Td) - Jva) to the line source function becomes impor- 
tant only in the transparent transition region (Fig. [12). In the 
real Sun locations with very bri ght H a, such as the moss and 
active-region heart in Fig. 13 of iRutt en (2007) may represent 
low-lying, denser transition regions that radiate Ha through re- 
combination paths. 



8. Conclusion 

The development of standard models of the solar atmosphere, 
masterminded by E.H. Avrett, represents a well-established pin- 
nacle of sophistication with respect to the application of NLTE 
spectrum formation theory with the inclusion of numerous spec- 
tral features. However, the assumption of hydrostatic and time- 
independent equilibria without magnetism remains a far cry 
from the actual solar chromosphere, which is pervaded by 
shocks and rapidly changing magnetic fine structure in even the 
quietest regions. 

Conversely, state-of-the-art magnetohydrodynamics simula- 
tions do a good job in emulating the small-scale magnetody- 
namism of the actual solar atmosphere, but they remain weak 
in properly treating non-equilibrium radiation. Implementation 
of the art of the ID modelers into the 3D time-dependent codes 
of the simulators presents a formidable challenge, but seems the 



|http : //www ■ staff . science . uu . nl/ - rutte IQ l/dot| 



most promisi ng venue to understan d the enigmatic solar chro- 
mosphere (cf . iLeenaarts et al]l2012lt . 
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